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Quantum technologies based on optical Gaussian states have proven very promising in terms of 
scalability. However, their use in quantum networking is hindered by the fact that Gaussian entan¬ 
glement cannot be distilled via Gaussian operations. We take advantage of hybrid optomechanical 
systems to address this problem, proposing a scheme to distill optical two-mode squeezed vacua 
via unsharp measurements. Here, one of the optical modes is injected into a single sided Fabry- 
Perot cavity and non-linearly coupled to a mechanical oscillator. Afterwards, the position of the 
oscillator is measured using pulsed optomechanics and homodyne detection. Our results show that 
this measurement can conditionally increase the initial entanglement under an optimal radiation- 
pressure interaction strength, which corresponds to an effective unsharp non-Gaussian measurement 
of the photon number inside the cavity. We show how the resulting entanglement distillation can 
be verified by using a standard quantum teleportation procedure. 


Introduction. — Recent experiments with quantum op¬ 
tics have demonstrated the generation of entanglement 
across hundreds of modes and partitions m, thus offer¬ 
ing unprecedented opportunities for quantum network¬ 
ing a- However, the states generated in these systems 
(Gaussian states of light) suffer from the drawback that 
entanglement distillation — a pivotal primitive for long 
distance quantum communication and networking |5 — 
is not readily available. This is due to the fact that 
the interactions naturally occurring in these systems are 
Gaussian and a “no-go theorem” prevents Gaussian op¬ 
erations to distill Gaussian entanglement |6]. In order to 
overcome this roadblock, purely optical methods involv¬ 
ing non-Gaussian operations have been suggested lam 
with the dominant scheme relying on photon subtrac¬ 
tion mi ng. The implementation of such schemes is 
currently topical but remains technologically challenging 
hms. We introduce here an alternative route based on 
hybrid opto-mechanical systems that naturally possess 
non-Gaussian radiation-pressure interactions. 

Quantum optomechanics is opening up new avenues 
for the manipulation of optical states P~9l [20] . The us¬ 
age of optomechanical systems for teleportation and es¬ 
tablishing Gaussian entangled states of distant systems 
have been studied (see, e.g. Refs. UMM)- How¬ 
ever the key quantum communication enabling proto¬ 
col of entanglement distillation has thus far been un¬ 
touched in opto-mechanics as the majority of the ap¬ 
plications considered a linearized (therefore Gaussian) 
interaction. On the other hand, the bare optomechan¬ 
ical radiation pressure interaction is non-Gaussian (tri- 
linear) [25]. Can this coupling be useful for quantum 
communication, namely, for entanglement distillation? 
The trilinear coupling, which enables mechanical super¬ 
positions [26H28] . has only recently started drawing seri¬ 
ous attention [29U37] as it is becoming physically signif¬ 



FIG. 1: (Color online) Distillation procedure for two-mode 
squeezed vacuum (TMSV). In a) we show the general scheme, 
where a source generates a TMSV, the first mode di interacts 
with a damped (/c) mechanical harmonic oscillator V ( x ) oc x 2 
(g is the scaled coupling strength). After this interaction, and 
conditioned on unsharp measurements of the position x of the 
mechanical oscillator with resolution S q , the initial TMSV en¬ 
tanglement between di and d <2 can be increased. In b) we show 
the explicit realization of the same procedure: here we sub¬ 
stituted the oscillator with a single sided Fabry-Perot cavity 
with a movable mirror, modeling the injection of the mode di 
into the cavity via BS of reflectivity r. The position of the 
mirror is to be measured by pulsed optomechanics S3]. 

icant in certain setups [38H43] . Here we show that the 
bare “tri-linear” optomechanical radiation pressure inter¬ 
action can enable the concentration of the entanglement 
of two mode squeezed vacua by local operations. In par¬ 
ticular, “snap-shot” position detections of a mechanical 
oscillator, whose technology has been developed recently 
[44d46] , serve as the alternative to photo-detection. Our 
proposal also demonstrates that sometimes the usage of 
weak (in the sense of “coarse-grained”/unsharp) mea¬ 
surements can be better for enacting a quantum protocol 
than their fine grained counterparts. 

System dynamics .— Let us commence by considering 
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two light-modes (with corresponding annihilation oper¬ 
ators di and <22 satisfying [dj,Oj\ = 1 for j = 1,2) in 
a two-mode squeezed vacuum (TMSV) |'0(O)) T msv = 
\J 1 — A 2 |n,n)i 2 , with A = tanhs and s being 

the squeezing parameter. One light beam (di) is coupled 
to a mechanical harmonic oscillator (ra), whereas mode 
a 2 propagates freely [a general scheme is illustrated in 
Fig. § a)]. As said, we focus our attention on a Fabry- 
Perot configuration [see Fig. §b)] where the mode d\ 
is injected into a cavity. Such injection of a propagating 
optical mode into a cavity is standard in LIGO @3, and 
is considered standard in the study of cavity-based quan¬ 
tum networks mmm- Obviously, the injection itself 
entails a decoherence of the optical held, which we will 
duly take into account. 


After mode di is injected into the cavity it starts inter¬ 
acting with the mechanical oscillator via an optomechan¬ 
ical Hamiltonian. In a frame rotating at the frequency 
co>i of mode di, this is given by Hint = b^b — gd\d\(b^ + 6), 
where g = go/oo m is the scaled coupling parameter, 
ujm is the angular frequency of the mechanical oscil¬ 
lator (b), go = x zp fUJi/L is the radiation-pressure in¬ 
teraction strength, L is the cavity length at equilib¬ 
rium, and x zp f is the zero-point fluctuation amplitude 
(we set h = 1) [26]. Given the recent possibilities 
of ground state cooling (231 SOI EMI] we will assume 
that the oscillator is initially in a coherent state | a) m . 
The evolution in the absence of any source of deco¬ 
herence can be solved straightforwardly 52. [53]. In 
this ideal case, the dynamics is characterized by a dis¬ 
placement of the mirror position, conditioned on the 
photon number n : A n | n) 1 \a) m A n exp(ig 2 n 2 (t — 
smt))exp(ignlm[arj]) \n) 1 \ae~ lt + gnrj ) . Here | n) is a 
photonic Fock state, g = 1 — e _2t , and t represents a 
scaled time, being the actual time multiplied by cj m . 


However, in realistic conditions the state will be af¬ 
fected by decoherence. In order to give a full analytic 
solution we assume that the cavity decay K> c /uj m « g 
(the resolved side-band regime already attained [54] , and 
a goal of several setups aiming to achieve ground state 
cooling). We solve the standard Markovian master equa¬ 
tion at zero temperature for the decoherence of the oscil¬ 
lator following the procedure in the appendix of Ref. |52j . 
In this case, the master equation reads as: dp(t)/dt = 
~i[H int ,p(t)] + (2bp(t)& - b ] bp{t) - p{t)&b)n/2, being k 
the mechanical energy damping rate. Another inevitable 
source of decoherence is the attenuation due to the injec¬ 
tion of the light beam into the cavity. To model this, we 
consider a beam splitter (BS) in front of the fixed cavity- 
mirror, such that one port of the latter is fed with mode 
d\ and the other with a vacuum held [55]. Under these 
sources of decoherence, the analytic solution for the full 


density matrix is 

00 
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where the 0 angle is related with the reflection coefficient 
of the BS as r = cos(0/2). The other terms are 
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Entanglement distillation .—To distill the initial TMSV 
we proceed to measure the quadrature position of the os¬ 
cillator [46] through an inefficient detector. This corre¬ 
sponds to the positive-operator valued measure (POVM) 
n(<?) = 1/V27r<5 g fZ exp (-(q - y) 2 /2<5 2 ) \y) (y\dy , 

where q = x \Jrnuj m /h is the dimensionless position of 
the oscillator (with actual position x), m is the oscilla¬ 
tor mass, and S q determines the measurement resolution 
156] . The state (unnormalized) after the measurement, 
conditioned to an outcome g, is given by 


minn.m 
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_ (q-x) 2 

where l nm = f™ 00 '4’<t> n ( K ,t)(x)'4’l rn ( K , t )(x)e 2S « dx , in 

which ^(g) = (g |£) is the position wave-function of 
an arbitrary coherent state |£). The probability density 
function (PDF) of the outcome g is 


p(9) 


|1~A 2 | 


y A 2 * exp 
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(q - V2Re [</>;(«, f)]) 2 
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(4) 

To quantify the entanglement we use the negativity [523- 
EU, defined as N(t) = 1/2 JV(|£*| — e*), where Si are the 
eigenvalues of the partial transposition of the normalized 
version of p(t) 12 of Eq. ©• 

A quick inspection of Eq. © reveals that a change 
in the initial amplitude from \a\e l( ^ a to \a'\e l( ^ a entails 
a rigid shift of the outcome probability p(q) by A q ~ 















3 


\a\ cos(</> a — t) — \a'\ cos(cj)' a — t). We verified numeri¬ 
cally that also the entanglement negativity is subjected 
to the same shift, which implies that a change in a can 
be accounted for by selecting the measurement outcome 
q accordingly. Given this, we set for the rest of this work 
the initial coherent state to a = 0. 

We now have the ingredients to assess the validity 
of the distillation procedure. For a fixed set of values 
(k = 0.01, <5^ ~ 0.11, r = 0.1, t = 7r, A = 0.3) we plot 
in the left y- axis of Fig. (|2]-a) the ratio of the negativity 
Nd/Nq (solid line) as a function of the outcome q of the 
measurement of the oscillator position, where Nd(Nq) 
stands for the distilled (initial) negativity. In the right 
y- axis, we show its corresponding PDF (dashed line) as a 
function of q. The success probability of the distillation 
protocol, namely the probability of obtaining Nd > We 
is given by the shaded region and is defined as: 

Pr (#,A) S = f p(q)dq, (5) 

J Nd>No 

In Fig. (|2|a) we illustrate three representative cases. For 
weak optomechanical coupling (g = 0.01), one achieves a 
large success probability though at the cost of an almost 
negligible increase in negativity Nd ~ Nq. For interme¬ 
diate coupling (g = 0.2) the negativity is significantly en¬ 
hanced, still retaining a high success probability. On the 
other hand, for large coupling (g = 1), not only Nd < Nq 
but also the success probability is considerably small. 

As a consequence, we see that an optimal region of 
the coupling value emerges, given that the entangle¬ 
ment distillation is predominantly achieved for inter¬ 
mediate radiation-pressure coupling. Similarly we can 
also see that entanglement distillation is achieved for 
intermediate values of the initial entanglement, imply¬ 
ing that distillation is optimal in the parameter region 
0.2 < {g, A} < 0.4, for which Pr(g,A) s > 0.15 and 
Nd > Nq. 

The reason for this behaviour can be intuitively un¬ 
derstood considering the structure of the TMSV state 
and its evolution under the distillation protocol. The 
states of the whole system (in absence of decoherence) 
before and after the optomechanical interaction are given 
b y l°) m En A ” K n >i ,2 and E„ \ n e i 9 2n2 ’ K |2 gn) m \n,n) 1>2 
respectively. The states \2gn) m become more and more 
distinguishable for larger g. As a consequence, the mea¬ 
surement of the oscillator position effectively becomes a 
sharp measurement of Fock state inside the cavity m 
that projects the two light beams into a factorized state 
|n, n)i ? 2 - This intuitively explains the failure of the dis¬ 
tillation protocol for large g. The failure for large A is 
instead due to the fact that the number of photon Fock 
states compatible with a specific outcome q is finite (for 
any non-zero g). For large enough A, this finite superpo¬ 
sition of a small set of Fock states |n, n)i,2 is not enough 
to exceed the entanglement of the initial TMSV. As said, 
we are neglecting here the cavity decay and further losses 



FIG. 2: (Color online) Upper panel a): We plot in the left 
?/-axis the ratio of the negativity Nd/No as a function of the 
oscillator’s position q (solid line), where Nd (No) stands for 
the distilled (initial) negativity. In the right y -axis we show 
the PDF as a function of q (dashed line). In the middle panel 
b) we illustrate the distillation success probability (Pr(g, A) s ) 
corresponding to the shaded region in the upper panel. Fi¬ 
nally, in the bottom panel c) we show the ratio of negativity as 
function of A and g for a specific oscillator’s position q — 1.5. 


in the extraction of the distilled state from the cavity. 
However, let us note that our results indicate that the dis¬ 
tillation protocol is robust against large injection losses 
(in Fig. [2] we considered a beam-splitter reflectivity of 
r = 0.1), which in turn suggests robustness against cav¬ 
ity and extraction losses as well. 


Quantum teleportation with the distilled state. —It is of 
course, crucial to suggest both a method to verify the suc¬ 
cessful distillation of entanglement, as well as an applica¬ 
tion of the distilled state. In the following we will show 
how the teleportation of an arbitrary coherent state \/3) 
by the distilled state can serve both purposes. Following 
the standard procedure m, we combine mode d\ with 
the coherent state to teleport into a balanced beam split¬ 
ter. Subsequently, we measure the position (momentum) 
quadrature of the transmitted (reflected) beam. The un¬ 
normalized state after the joint {x,p}-measurements cor- 
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FIG. 3: The figure shows the ratio of the average fidelity 
({Ed(x,p)) p / {Eo(x,p)) p) as a function of the amplitude (/3) 
and phase (c/)p) of the coherent state to teleport. In this figure 
we used the parameters a — 2e Z7F//4 , q = 0, A = 0.3, g — 0.2. 
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where V(n — k,j, n',j') = ( n j fe ) (^e^C 2 / n,) H j+j , (x) x 
iJ n _fc +n /_j_j/(p), in which the Hermite polynomial 
of degree n. The normalization of Eq. ([6| gives us the 
probability of the joint measurement p(x,p) = Trp(t) 2 - 
The success of the teleportation can be quantified by the 
fidelity which in our case reads as f(x,p) = (/3\p2 |/3), 
being p 2 the normalized state appearing in Eq. © dis¬ 
placed by D(y/2(x + ip)). In particular, we consider the 
fidelity averaged over all possible outcomes: 


/ +oo 

f(x,p)pp(x,p)pdxdp. (7) 

-OO 

In Fig. ([ 3 ]) , with parameters a = 2e Z7r / 4 and q = 0, 
we clearly distinguish many instances in which the fi¬ 
delity with the distilled state is greater then the one with 
the original state, z.e., > {Fo{xiP))p [the 

subindex in Tb (T§) corresponds to the distilled (initial) 
state]. Thus, for a wide range of /3, the distilled state can 
be exploited to improve teleportation fidelity. 

Experimental feasibility. — The generation of the 
TMSV state is routine in quantum optics laboratories 
where 0 < A < 0.5 can be fully realized. On the other 
hand, our protocol requires the initialization of the me¬ 
chanical oscillator in its ground state or in a coherent 
state. The cooling of a mechanical oscillator to its quan¬ 


tum ground state has been recently achieved [40l|5T]. Fi¬ 
nally, typical values in optomechanics show a wide range 
for the oscillator’s damping rate k [39, 4TJ [62 compati¬ 
ble with our assumptions. Moreover, our optimal g value 
of ~ 0.2 for the optomechanical interaction can now be 
achieved in several setups [38H43] . 

The key stage of this work consists in measuring the 
oscillator. After the pulse di interacts with the oscilla¬ 
tor under sufficiently weak coupling (e.g., using g ~ 0.2, 
cj m /27r = 500 kHz, and go ~ 86 kHz, as in |56|), a second 
auxiliary pulse with a duration much smaller than the 
mechanical period is then injected into the cavity. The 
optical phase of the emerging field (correlated with the 
mechanical position) is then measured via balanced ho¬ 
modyne detection [46]. Moreover, ideally, we want our 
optomechanical coupling to cease after the measurement 
as we do not want entanglement of the distilled quan¬ 
tum state with mechanics while it is being used for some 
other protocol. This can be achieved, for example, by 
using trapped mechanical objects as oscillators interact¬ 
ing with a cavity field: an optically trapped object can 
be pushed away by a strong pulse [63]; a charged object 
held in a Paul trap can be suddenly shifted relative to 
the cavity field by suddenly shifting the trap centre [64, . 

Concluding remarks. — We have presented a first appli¬ 
cation of the bare radiation-pressure coupling in a prac¬ 
tical quantum communication protocol, showing how op¬ 
tomechanics can fill a crucial gap in enabling long dis¬ 
tance quantum networks using Gaussian states of light. 
Our proposal uses an indirect measurement of the photon 
number of the field inside a cavity through the position 
measurement of a mechanical element coupled to it. For 
an optimal strength of the coupling, the photon number 
is measured weakly or “unsharply” and this results in en¬ 
tanglement distillation conditioned to the position out¬ 
come. For a vacuum state inside the cavity, the position 
meter does not move, corresponding to a failed outcome 
of distillation. Thus our procedure has a degree of simi¬ 
larity with the known purely optical procedure of photon 
subtraction [MU, where also the vacuum component is 
filtered, and we have comparable figures of merit. On the 
other hand, the position measurements of a mechanical 
oscillator can be highly precise, especially as we are not 
concerned about back-reaction as the measurement is at 
the end of our protocol (the oscillator can be reinitial¬ 
ized before performing distillation again). Optimizations 
with multiple modes in the cavity can be attempted. 

The state obtained through our protocol is non- 
Gaussian, and thus it can serve as the first step of Gaus- 
sification jgiss] — which can enhance its entanglement 
further and act on multiple copies — or, more in general, 
for quantum computation purposes [66] . Moreover, the 
procedure here outlined could be useful also in a quantum 
repeater scenario for long distance communication, con¬ 
sidering that there no further extraction of the distilled 
state from the optical cavity is needed. 
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